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We report results of the Dalitz analysis of the three-body charmless B + — > K + tt + tt~ and B+ — * 
K+K+K~ decays based on a 140 fb™ 1 data sample collected with the Belle detector. Measurements 
of branching fractions for quasi-two-body decays to scalar-pseudoscalar states: B+ — > fo(980)K + , 
B+ -» if*(1430)°7r+, and to vector-pseudoscalar states: B+ -» if*(892)°7r+, B+ -> p°K+, B+ -* 
tj>K+ are presented. Upper limits on decays to some pseudoscalar-tensor final states are reported. We 
also report the new measurement of the _B+ — > x c oK+ branching fraction in two XcO decay channels: 
XcO —* 7r+7r~ and \c0 — > K + K~ . 



1 Introduction 

Studies of S meson decays to three-body 
charmless hadronic hnal states are a natu- 
ral extension of studies of decays to two- 
body charmless final states. Some of the final 
states considered so far as two-body (for ex- 
ample, pTt, K*tt, etc.) proceed via quasi-two- 
body processes involving a wide resonance 
state that immediately decays in the sim- 
plest case to two particles, thereby produc- 
ing a three-body final state. B meson decays 
to three-body charmless hadronic final states 
may provide new possibilities for CP viola- 
tion searches. 

Observation of B meson decays to var- 
ious three-body charmless hadronic final 
states has already been reported by the 
Belle 1 ' 2 ' 3 , CLEO 4 and BaBar 5 Collabo- 
rations. First results on the distribution of 
signal events over the Dalitz plot in the three- 
body B+ -> K+tt+tt- and B+ -> K+K+R- 
decays are described in Ref. 1 . With a data 
sample of 29.1 fb~ x a simplified analysis tech- 
nique was used because of lack of statistics. 
Using a similar technique, the BaBar collabo- 
ration has reported results of their analysis of 
the Dalitz plot for the decay B + — > K+tt+tt^ 
with a 56.4 fb" 1 data sample 6 . With the 
large data sample that is now available, we 
can perform a full amplitude analysis. The 
analysis described in this paper is based on a 



140 fb^ 1 data sample containing 152 million 
BB pairs, collected with the Belle detector 
operating at the KEKB asymmetric-energy 
e + e~ collider. 



2 Amplitude Analysis 

Analysis of two-body mass spectra shows 
that a significant fraction of the signals ob- 
served in B + — > K+tt+tt^ and B + — > 
K + K + K~ decays can be assigned to quasi- 
two-body intermediate states. These reso- 
nances will cause a non-uniform distribution 
of events in phase space that can be analyzed 
using the technique pioneered by Dalitz. 
Multiple resonances that occur nearby in 
phase space will interfere and provide an op- 
portunity to measure both the amplitudes 
and relative phases of the intermediate states. 
This in turn allows us to deduce their relative 
branching fractions. Details of the event se- 
lection and amplitude analysis could be found 
in the Belle contributed paper 7 . Here we 
can present the main results only. The ex- 
amples of the two-body invariant mass dis- 
tributions and their description by a fit of 
B+ -> K+tt+tt- and B+ -> K+K+K~ de- 
cays are presented in Fig. 1. Results on the 
branching fractions for quasi-two-body de- 
cays are summarized in Table 1. 
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Figure 1. Results of the fit to _R"+7r+7r and K+K^K events. Points with error bars are data, the open 
histograms arc the fit result and hatched histograms are the background components. Insets show: (a) the 
K*(892) - JC*(1430) mass region in 20 McV/c 2 bins; (b) the \c0 mass region in 25 McV/c 2 bins; (d) the 
0(1020) mass region in 2 MeV/c 2 bins. 



Table 1. Summary of branching fraction results. The first quoted error is statistical, the second is systematic 
and the third is the model error. The charmless total fractions in this table exclude the \c0 contribution. The 
value given in brackets for the Kq(1430)7t+ and XcoK + channels corresponds to the second solution (see 7 
for details). 



Mode 


B(B+ -» Rh+) 
xB{R->h+h~) x 10 6 


B(B+ -» Rh+) x 10 6 


K + ir + ir~ charmless total 




46.6 ±2.1 ±4.3 



K*(892)°ir+, K*(892)° -> K+n- 
X *(1430)7r+, #5(1430) K+n~ 

^*(1410)tt+, A"*(1410) K+n- 
K*(1680)ir+, X*(1680) K+Tr~ 
i^ 2 *(1430)7r+, #|(1430) K + tt~ 
p°(770)K+, p°(770) tt+tt- 
/o(980)#+, / (980) -» tt+tt- 
/ 2 (1270)A' + ) / 2 (1270) - tt+tt" 
Non-resonant 



6.55 ±0.60 ±0.60^;^ 

27.9 ± 1.8 ± 2.6±|;5 
(5.12 ± 1.36 ± 0.49±o.li) 

< 2.0 

< 3.1 

< 2.3 

4.78 ± 0.75 ± 0.44±g$ 
7.55 ± 1.24 ± 0.69±o;|l 

< 1.3 



9.83 ± 0.90 ± 0.90±g;|^ 
45.0 ± 2.9 ± 6.2± 1 |;^ 
(8.26 ± 2.20 ± 1.19±g;g|) 



4.78 ± 0.75 ± 0.44±g;|^ 



17.3 ± 1.7 ± 1.6±\ 7 8 



K + K + K charmless total 
4>K+, 4>-> K+K- 
0(1680)^+, 0(1680) -> K+K- 
/o(980)#+, / (980) -> K+K- 
ft(1525)K+, $(1525) -» 
02(1320)^+, a 2 (1320) -> #+#" 
Non-rcsonant 



4.72 ± 0.45 ± 0.35±g;H 

< 0.8 

< 2.9 

< 2.1 

< 1.1 



30.6 ±1.2 ±2.3 
9.60 ± 0.92 ± 0.71±g;I| 



1.37 ± 0.28 ±0.12+^| 
0.86 ± 0.26 ± 0.06^ 05 
(2.58 ±0.43 ±0.19^05) 



24.0 ± 1.5 ± 1.8±^ 



Xc0K + , XcO ~> 7T + 7T 

XcoK+, XcO -> ^+^" 
X c o# + combined 



196 ± 35 ± 33laf 



3 Discussion 

With a 140 fb _1 data sample collected with 
the Belle detector, we have performed the 



first amplitude analysis a of B meson dc- 

a At this Conference Babar presented the Dalits plot 
analysis for K + n + ir~ 8 . 
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cays to the three-body charmless K + Tr + n~ 
and K + K + K~ final states. Such analysis of- 
ten suffers from uncertainties related to the 
non-unique description of the decay ampli- 
tude. In our case the uncertainty comes 
mainly from the parametrization of the non- 
resonant amplitude. It is worth noting that 
fractions of the non-resonant decay in both 
B+ K+tt+tt- and B+ K+K+K~ de- 
cays are comparable in size and comprise a 
significant fraction of the total three-body 
signal, which may indicate the common na- 
ture of the amplitudes. 

Despite the large model uncertainty dis- 
cussed above, there is a set of quasi-two- 
body signals whose branching fractions can 
be measured with a relatively small model 
error. In particular, clear signals are ob- 
served in the B+ if*(892)°7r+, B+ -> 
p°(770)K+, B+ f (980)K+ and B+ -> 
(f)K + decay channels. The model uncertainty 
for these channels is small due to the nar- 
row width of the resonances and in vector- 
pseudoscalar decays due to the clear signa- 
ture of the vector meson polarization. 

The branching fraction value measured 
for the decay B+ -> Jf*(892)°7r+ is signif- 
icantly lower than that reported earlier 1 ^ 6 . 
The simplified technique used for Dalitz anal- 
ysis of the B + — > K + tt + tt~ decay described 
in 1,6 has no sensitivity to the relative phases 
between different resonances, resulting in 
a large model error. The full amplitude 
analysis presented in this paper consistently 
treats effects of interference between differ- 
ent quasi-two-body amplitudes thus reducing 
the model error. The analysis suggests the 
presence of a large non-_ft"*(892)°7r + (presum- 
ably non-resonant) amplitude in the mass re- 
gion of the K* (892)° that absorbs a signif- 
icant fraction of the B signal. The B + — > 
K*(892)°ir + branching fraction measured in 
our analysis is in better agreement with theo- 
retical predictions based on the QCD factor- 
ization approach. 

The decay mode B+ / (980)^+ 



is the first observed example of a B de- 
cay to a charmless scalar-pseudoscalar fi- 
nal state. The mass M(/ o (980)) = 976 ± 
4± 2 MeV/c 2 and width r(/ (980)) = 61 ± 
9^g MeV/c 2 obtained from the fit arc in 
agreement with previous measurements. The 
sensitivity to the B + — ► fo(980)K + decay in 
the K + K + K~ final state is greatly reduced 
by the large B + — ► (f>K + signal and by the 
scalar non-resonant amplitude. No statisti- 
cally significant contribution from this chan- 
nel to the K + K + K~ three-body final state 
is observed, thus only a 90% confidence level 
upper limit for the product of the correspond- 
ing branching fractions is reported. 

We report the first observation of the de- 
cay B+ p°(770)K + . This is one of the 
channels where large direct CP violation is 
expected. 

Due to the very narrow width of the <j> 
meson, the branching fraction for the decay 
B + — > <pK + is determined with a small model 
uncertainty. 

A clear signal is also observed for the 
decay B + — > xmK + in both xm —* tt + tt~ 
and XcO ~ * K + K~ channels. Although quite 
significant statistically, the B + — > XcoK + 
signal constitutes only a small fraction of 
the total three-body signal and thus suffers 
from a large model error, especially in the 
K + K + K~ final state. For this decay mode, 
the charmless non-resonant amplitude in the 
XcO mass region is enhanced compared to the 
K + tt + tt~ final state due to the interference 
caused by the presence of the two identical 
kaons. 

For other quasi-two-body channels the 
interpretation of fit results is less certain. Al- 
though the B+ -> if *(1430)7r+ signal is ob- 
served with a high statistical significance, its 
branching fraction is determined with a large 
model error. Two solutions with significantly 
different fractions of the B+ X£(1430)tt + 
signal but similar likelihood values arc ob- 
tained from the fit to K + tt + tt~ events. A 
study with MC simulation confirms the pres- 
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ence of the second solution. This may indi- 
cate that in order to choose a unique solution 
additional external information is required. 
In this sense, the most useful piece of infor- 
mation seems to be the phenomenological es- 
timation of the B+ -> ^(1430)tt + branch- 
ing fraction. The analysis of B meson decays 
to scalar-pseudoscalar final states described 
in Ref. 9 suggests that the branching fraction 
for the B + — > Kq (1430)7r+ decay can be as 
large as 40 x 1CP 6 . Unfortunately, the pre- 
dicted value suffers from a large systematic 
error that is mainly due to uncertainty in the 
-?£Tq(1430) decay constant fK*(U30)- Differ- 
ent techniques used to estimate (1430) 9 ' 10 
give significantly different results. Further 
improvement in this field would be useful. 

We also check possible contributions 
from B decays to various pseudoscalar-tensor 
(PT) states. In the factorization approxima- 
tion, charmless B decays to PT final states 
are expected to occur at the level of 10 -7 or 
less. We find no statistically significant sig- 
nal in any of these channels. As a result, we 
set 90% confidence level upper limits for their 
branching fractions. 

We cannot identify unambiguously the 
broad structures observed in the M(ir + ir~) ~ 
1.3 GeV/c 2 mass region in the K + ir + ir~ fi- 
nal state denoted as /x(1300) in our analy- 
sis and at M{K + K~) ~ 1.5 GcV/c 2 in the 
K+K+K- final state denoted as / x (1500). 
If approximated by a single resonant state, 
/x(1300) is equally well described by a scalar 
or vector amplitude. Analysis with higher 
statistics might allow a more definite conclu- 
sion. The best description of the /x(1500) 
is achieved with a scalar amplitude with 
mass and width from the fit consistent with 
/o(1500) states. 

Results of the B+ K+K+R- Dalitz 
analysis can be also useful in connection with 
the measurement of CP violation in B° — > 
KgK + K~ decay reported recently by the 
Belle 11 and BaBar 12 collaborations. An 
isospin analysis of B decays to three-kaon 



final states suggests the dominance of the 
CP-even component in the B° — ► KgK + K~ 
decay (after the B° — ► <f>Kg signal is ex- 
cluded). This conclusion can be checked in- 
dependently by an amplitude analysis of the 
KgK + K~ final state, where the fraction of 
CP-odd states can be obtained as a fraction 
of states with odd orbital momenta. Unfor- 
tunately, such an analysis is not feasible with 
the current experimental data set. Neverthe- 
less, the fact that we do not observe any vec- 
tor amplitude other than B + — > <fiK + in the 
decay B + —> K + K + K~ confirms the conclu- 
sion. 
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